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The Advanced Receiver currently being developed uses a Costas digital loop 
to demodulate the subcarrier. Previous analyses of lock detector algorithms for 
Costas loops have ignored the effects of the inherent correlation between samples 
of the phase-error process. Accounting for this correlation is necessary to achieve 
the desired lock-detection probability for a given false-alarm rate. In this article, 
both analysis and simulations are used to quantify the effects of phase correlation 
on lock detection for the square-law and the absolute-value type detectors. Resui s 
are obtained which depict the lock-detection probability as a function of loop signal- 
to- noise ratio for a given false-alarm rate. The mathematical model and computer 
simulation show that the square-law detector experiences less degradation due to 
phase jitter than the absolute-value detector and that the degradation in detector 
signal-to-noise ratio is more pronounced for square-wave than for sine-wave signa s. 


I. Introduction 

Costas loops are being used extensively in modern co- 
herent communication systems to track both sub carriers 
and suppressed carriers. In many applications, residual 
carriers are being replaced by suppressed carriers as the 
latter dedicates the total transmitted power to both car- 
rier tracking and symbol detection simultaneously. This 
has a clear advantage over residual carrier tracking, which 
requires a fraction of the total transmitted power delegated 
solely to that purpose, and hence, reduces the available 
power that can be used for symbol detection. The dis- 
advantages of Costas loops are that they require symbol 
synchronization and suffer from an additional loss factor 


typically referred to as squaring loss, which is highly depen- 
dent on symbol energy-to-noise ratio. Squaring loss is the 
result of forming the product of the inphase and quadra- 
ture signals to wipe out the data modulation, in order to 
obtain a feedback error signal that is only a function of 
the instantaneous phase error [1]. Another disadvantage of 
suppressed carrier tracking is the issue of false lock, which 
occurs either as a result of accumulated delay [2] or during 
acquisition with frequency uncertainty greater than one- 
half the symbol rate [3, 4]. The latter can be detected 
through a false-lock indicator, as described in [3]. 

Lock detection is an important part of a tracking 
loop’s operation and monitoring, as it provides an insight 
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into the tracking loop’s behavior in real time. Lock de- 
tection basically serves as a binary indicator of whether 
the loop is tracking the received signal or not, and during 
loop start-up it also indicates whether or not the loop has 
acquired the phase of the signal. There are mainly two 
kinds of lock detectors for Costas loops, the I 2 — Q 2 , or 
square-law detector, and the | I \ — \ Q |, or absolute- value 
detector [5]. Both have been analyzed in the past at high 
loop signal-to-noise ratio (SNR), which basically assumes 
zero phase-tracking error. At low loop SNR, the assump- 
tion of zero phase jitter becomes inadequate and results in 
system operating parameters that are different from their 
design counterparts. Thus, a new model is required which 
has to account for the phase jitter and the correlation be- 
tween samples of the phase-error process. The latter is es- 
sential in an accurate performance prediction analysis, in 
order to operate the system at the desired lock-detection 
probability for a given false-alarm rate. The analysis of 
the detectors, including the phase correlation and assum- 
ing either a sinusoidal or a square wave signal, is presented 
in Section II. In the case of a square- wave, the analysis is 
general and includes any windowing operation as described 
in [6]. The discussion of the results and some simulated 
data are shown in Section III, followed by the conclusion 
in Section IV. 


II. Lock-Detection Analysis 

Suppressed carrier tracking for binary-phase shift 
keyed (BPSK) signals can be accomplished using a squar- 
ing loop or a Costas loop [7]. The squaring loop relies on 
wiping out the data by a squaring operation and tracking 
the resulting residual double-frequency component with a 
classical phase-locked loop. The squaring loop does not 
require symbol timing, but results in an additional noise 
term which becomes dominant at low SNR. On the other 
hand, the Costas loop implements a phase discriminator 
by forming the product of the inphase and quadrature sig- 
nals. That too results in some degradation, commonly 
referred to as the squaring loss. Depending on the de- 
sign, both loops can be implemented with identical per- 
formances for all practical purposes. The Costas loop has 
various derivatives, each approximating the maximum a 
posteriori (MAP) estimator at different SNRs [8]. For ex- 
ample, a hard-limiter can be included in the inphase arm 
to estimate the current symbol, and that results in less 
squaring loss at high SNRs. 


dump arm filters are digital accumulators. The IQ loop 
and square-law lock detector are depicted in Fig. 1 for the 
square-wave case, with the optional windowing operation 
on the quadrature channel. The analysis that follows will 
be applicable for both sinusoidal and square waves. The 
received waveform is digitized to produce the samples rj, 
which are subsequently digitally mixed with the in phase 
and quadrature references. The outputs of the mixers, 
running at the sampling rate, are accumulated to detect 
the received symbols. It is assumed that there are L sam- 
ples per symbol and that perfect symbol synchronization 
has been achieved. The accumulator outputs, now at the 
symbol rate, are multiplied together to wipe out the data 
and again accumulated to reduce the processing rate even 
further. The output, running at a new rate (referred to 
as the loop update rate), is the input to the digital loop 
filter which provides a frequency estimate to adjust the 
phase of the numerically controlled oscillator (NCO). The 
lock detector processes the arm accumulator outputs at 
the symbol rate, accumulates the result over M symbols, 
and provides a binary decision on the loop status. 

In the I and Q branches of Fig. 1, the signals accu- 
mulated over L samples during the k th symbol interval are 
given by 


xj k — d k L\J Pd Wk + Tijk fc — 

and (1) 

%Qk = dkL\ / Pd v k + TiQk k = 1, . . . , M 

where 

2 

Ufc = ~<f> k, I <f>k \< w 

7T 

I 4k |< *W/2 ( 2 ) 

•kW/2 <| <f>k |< tt(1 — T^/2) 
l 2sgn(<£it) - u* 7r(l - W/2) <| 4k |< it 


wt = (l- I |), 


Vk = < 


sgn (4 k)W 


for a square- wave sub carrier and 


u k - COS <f) k 


This article is concerned with the lock detection for 
the all-digital IQ loop, which is also a derivative of the 
Costas loop with integrate-and-dump arm filters. All- 
digital refers to the fact that the input waveform to the 
loop is a sequence of samples and that the integrate-and- 


Vk 


( sin^jb | <f) k |< itW/2 
\ sm(<j> k lV) nW/2 <\ <j> k |< tt( 1 - W/2) 


( 3 ) 


t sin <f>k tt(1 - W/2) <| 4 k |< x 
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and 


for a sine- wave subcarrier, where Pd is the average data 
power, djfc is the data value of the kth binary symbol (±1 
equally likely), (j>k is the subcarrier phase-estimation error 
(radians) at time Jb, W is the width of the window in the 
Q channel (i.e. , the fraction of cycle of the reference signal 
which has nonzero value W < 1 [W = 1 means no window 
is used]), and TiikifiQk are zero- mean white Gaussian noise 
samples. From Eq. (1)? the mean values and the variances 
of xn and XQk conditioned on 4>k and dk are given by 

fJ'Ik ~ dkL\fPo cr? = L a n 

and ( 4 ) 

VQk = dkLs/Po Vk &q = WLcr 2 

where <j\ = N 0 B n is the noise variance of a received sam- 
ple, N 0 is the one-sided noise spectral density, and B n is 
the Nyquist bandwidth. The equations that follow are ap- 
plicable to both sine and square waves. 

A. Square-Law Detector 

The first algorithm considered detects the in-lock state 
by producing a signal that is proportional to the cosine 
of the phase error (in the case of a sinusoidal wave) and 
averaging it over several symbols before comparing it to a 
threshold r. Referring to Fig. 1, 

M 

y* < r ’ wliere yk = _ x h k 

fc=l 

An estimation of the performance of this lock detector re- 
quires the first and second moments of and con- 
ditioned on <pk^ which are readily obtainable: 


x n “ P/fc + a i (^ a ) 

^ = + 4 ( 6b ) 

xj k = n) k + SpM + 3 crj (7a) 

X Qk = ^Qk + $l*Qk a Q + 3^0 ( 7l3 ) 


var(xg fc ) = 4nQ k o-Q + 2a q (8b) 

The mean value of y k is obtained from Eqs. (4) and (6) in 
Eq. (5), namely 

liy k = L 2 P D {wl-vl) + L**{ l-W) 

The variance of y k , (Ty k , will be the sum of variances of x j k 
and XQ k , and is obtained by using Eq. (4) in Eq. (8), i.e., 

< = + Wv D + 1 + 

The lock-detector signal z is the accumulation of Myk sam- 
ples, which are highly correlated due to the phase-error 
samples. The mean value and variance of z, \i z , and 
are derived in the Appendix, where it is shown that 

p z = m(l 2 P d < 1+ Lal(l-wj) ( 9 ) 

and 

a 2 = M 2 L 4 Pj)(g - d 2 ) + 4 ML 3 P D a 2 n (f + hW) 

+ 2ML 2 a 4 (l + W 2 ) ( 10 ) 

The parameters d, f, h, and g depend on the waveform 
type and are given by (b = Wtt/2) 

— * (!)'- 

Q / A\ 2 1 M_1 

0 = i - - 1 + (-) jp 12 c ( k ) d ( k ) 

v ' k - 0 


where y Ik , HQk, 0j, and a 2 Q are given by Eq. (4). Using 

Eqs. (6) and (7), the variances of X/*, and Xg*. become where 

respectively 


var(z/ fc ) = ($ a ) 
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for k = 0 


and 


for a square wave and 


<t> 2 ~ a\ erf 



d{k)={ 


0.5(l + exp(-2<r£)) 
exp(— 4^) cosh (2<T^C(T k )j 


fc = 1,2, . . . , Af - 1 



for a square- wave sub carrier, and 
d = exp(— 2cr|) 


for a sine wave; C(rk) is the correlation function of the 
phase-error process in the tracking loop, which is assumed 
to be of the form given by Eq. (A-9) [9], and the second- 
order joint probability density function of the correlated 
phase process p{<f>u <f>2) 7 *) is assumed as in Eq. (A-10). 
When the loop is in-lock and assuming high loop SNR, 
<p k — ► 0 for all k . Hence 1 and v\ — ► 0, and the above 

mean value and variance of the detector’s signal simplify 
to 


/ = 0.5 (1 +exp(— 2<r£)) 


fi z = m(l 2 p d + L*l( 1 - wo) 


h = 0.5(l - exp(— 2cr|)) 


a 2 = AML 2 at 


LP d 1 + W 2 


+ 


9 = 


1 

M 2 


M - 1 

E 


fc = 0 


which are true for both square-wave and sine-wave sig- 
nals. Note that in the above equation the following rela- 
tion holds: 


2 _ f Nq B sc 

V Pd 


Wi 1 + 


2E S /Nq 


( 12 ) 


LPd 2LL 


for a sine- wave sub carrier. Note that Eqs. (11) and (12) 
specify the variance of the phase jitter, assuming the linear 
loop model. For example, at 15 dB of nominal loop SNR, 
the actual variance can be about 1 dB larger. Nominal 
loop SNR p is defined as 1/trJ, where is obtained from 
the linear model; B sc is the one-sided noise bandwidth of 
the Costas loop, and E s /Nq is the symbol energy- to-noise 
ratio. The constants c(k ) and d(k ) are defined as follows: 

f M for k — 0 

c(Jfc) = l 

v y [2{M-k) for k — 1,2,...,M — 1 
for both waveforms, but 


because L = 2 T s B n (Nyquist sampling), where T s is the 
symbol time, <r\ = N 0 B n , and T S P D = E S} the energy per 
symbol. 


B. Absolute-Value Detector 

For the absolute- value detector, the squaring opera- 
tion in Fig. 1 is replaced by an absolute- value operation, 
and the algorithm defining the new lock detector becomes 


M 

5^ Vk \ T y where y k = \ x Ik | - | x Qk \ (13) 

fc= 1 


b b 

d(k)= j j | <f>l || <f>2 | 

-b -b 


In order to estimate the performance of this lock detector, 
the first and second moments of | xj k | and | XQ k | are 
needed, again assuming a white Gaussian noise process at 
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the phase-locked loop input. These moments, conditioned 
on <j>k, become 

1 T A 

I zj* I = r n 



phase process in the tracking loop. This is carried out in 
the Appendix and gives 

fi z - M fly 





a 2 = m(l 2 P d (w* + u 2 ) + La 2 n { 1 + W)) 


+ + r Q iT Q^ - M 2 ( ri 2 + r< 5 2 ) 

all ij 

(17) 


The second moments of | x/jc | and | zq* | are identical 
to the second moments of x/ k and XQk, and are given by 
Eq. (6). The mean value of y k follows from Eqs. (13) and 
(14): 

A 

Hyk 7 r Ik -- r Qk 



(15) 


and the variance of y k is 

< 7 2 k = L 2 P D (ul + V 2 ) + La 2 (l + 17) - (m 2 + TQk 2 ) 

The lock detector’s signal z is again obtained by adding 
My * samples. The mean and variance of z are found by 
averaging the first two moments of z over the correlated 


The bar over the product terms denotes expectation over 
the joint probability density function of <pi , <p j , assumed 
to be of the form given by Eq. (A-10), and r Ifc , rq k are 
defined by Eq. (14). Because no closed-form solutions for 
the above averaging operations are known, the averaging 
was done numerically. When the loop is in-lock and at 
high loop SNR, <pk — ► 0 for all k. Hence, u> k ♦ 1 and 
— > 0 } and the above mean value and variance of the 
lock detector simplify to 



<r 2 = M L 2 P d + L<r 2 n (l + W) 
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C. Probability of Detection and of 
False Indication 

During subcarrier detection, each z sample is com- 
pared with a predefined threshold r, and the lock detector 
decides that the loop is in-lock when z > r . It is possible 
that even when no signal is present, z will occasionally be 
larger than r. In this case, the lock detector will mistak- 
enly declare an in- lock condition. The probability of false 
indication is 


Pf 


1 

V^To 



{z -Hto) 2 \ 
/ 


dz 



r- fizo \ 

rjoj 


( 20 ) 


where \i z o and <t^ 0 are the mean and variance of the lock- 
detector signal in the out-of-lock state, and erfc(x) is the 
complementary error function (erfc(x) = 1 — erf(x), where 
erf(x) is the error function defined in the Appendix). For 
the square-law detector, fi z o and cr] 0 are obtained from 
Eqs. (9) and (10) by making P D = 0 (or, equivalently, 
assuming that UT^ and vj in Eq. 4 are zero), namely 


VzQ = ML<t£(1 — W) 

<j\ 0 — 2ML 2 <t*(1 + W 2 ) 

whereas for the absolute- value detector, Eqs. (18) and (19) 
result in 



Given a desired probability of false indication P j , the 
threshold r is obtained by solving Eq. (20) and setting 
it equal to 


t — \J 2<r^ 0 erfc 1 (2P/) + p z o (22) 

where erfc -1 () is the inverse complementary error func- 
tion. When the loop is in-lock, it can be argued via the 
central limit theorem that the random variable z is ap- 
proximately Gaussian, with mean and variance obtained 
earlier for either the square-law or the absolute-value de- 


tector. For either detector, the probability of detection 
is 





2*1 J 


dz 





where p z and <r^ are given by Eqs. (9) and (10) or by 
Eqs. (16) and (17). Defining the detector’s SNR as 


SNR; 



(23) 


then, for p z o = 0 ( W — 1), the probability of detection in 
terms of SNR; can be expressed as 


Pd = 2 erfc 


— erfc" 1 (2 P t )~ 

O’z 


SNR; 


The above equation shows the dependence of the proba- 
bility of detection on the detector’s SNR. Phase jitter in 
the tracking loop degrades the detector’s SNR by a factor 
D: 


D = SNR; /SNR; (id ea r) (24) 

where SNR; (,•*„,) is the detector SNR, assuming infinite 
loop SNR, i.e., no phase jitter. SNR; ( i<jea( ) is computed 
from Eq. (23) using the high-SNR expressions in Eqs. (18) 
and (19) for p z and a 2 z . For a given W, M, loop SNR p 
(p = l/tf^ where is given by Eq. 11 or 12 repectively), 
and Pf, the detector’s SNR must be increased approxi- 
mately by the factor 1/D in order to achieve a desired 
probability of detection. 


III. Discussion and Numerical Results 

Computer simulation was performed in order to check 
the predictions of the analysis. Figure 2 depicts the prob- 
ability of lock detection versus symbol energy-to-noise ra- 
tio E s /Nq for both sine-wave and square-wave signals, as- 
suming ideal conditions, i.e., no phase jitter in the track- 
ing loop. The square-law detector performs slightly better 
than the absolute-value detector for a given symbol SNR. 

The degradation in detection probability for a finite- 
loop SNR is shown in Figs. 3(a) and 3(b) for the square-law 
and absolute-value detectors, respectively. The threshold 
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r wets set to achieve probabilities of false detection Pj of 
1CT 1 and 10“ 4 , and detector SNR was set to achieve nomi- 
nal probabilities of detection P d of 0.99 and 0.90. Nominal 
probability refers to the case of no phase jitter in the loop. 
It is clear that sine waves produce less degradation than 
square waves, and this is true for both detection schemes. 

The performance of both detectors is compared in 
Fig. 4 for square-wave signals only, since the difference 
in performance is almost negligible for sine-wave signals. 
The performance with respect to detector SNR is shown 
in Fig. 5 for a 15-dB loop SNR. The improvement in de- 
tection probability due to windowing is clear for both de- 
tectors and can result in several decibels. Finally, Fig. 6 
depicts both theoretical and simulation points of detec- 
tor SNR degradation D (defined in Eq. 24) versus loop 
SNR. The degradation in SNR is slightly larger for the 
absolute-value detector than for the square-law detector 
when tracking a square wave, but less when tracking a 
sine wave, and it can be as large as 3 dB depending on the 
operating parameters. 

The results are summarized in Figs. 7 and 8. The 
detector SNR as a function of E s /N 0 is shown in Fig. 7 
for both infinite and 15-dB loop SNR, and for W = 1.0 
and 0.25. When W = 1, a good rule of thumb is that 
the detector SNR varies linearly with E 9 /Nq , with slope 
equal to 2/3 on a decibel scale. For different values of M , 
the curve will be scaled vertically in a linear fashion. Fig- 
ure 8 depicts the detection probability for the square-law 
and absolute-value detectors respectively, as a function of 
E s /Nq for both infinite and 15-dB loop SNR. The con- 
clusion from Fig. 8 is that when operating at low loop 
SNR (i.e., 15 dB), an extra 1.5-dB increase in E s /Nq or a 
comparable increase in M will achieve the detection prob- 
ability which was designed for assuming infinite loop SNR. 

For design purposes, Fig. 9 can be useful since it 
depicts both the detection probability and the required 


threshold as a function of M for both detectors. As a 
design example, suppose that the absolute-value detector 
is required to operate at Pj = 10 -4 and Pd = 0.99, and 
that the signal is a square wave with symbol rate r s = 80 
symbols per second, E s /Nq = 0.0 dB, and loop SNR = 
15 dB with a quarter window ( W = 0.25). Figure 9 indi- 
cates that at least 90 detector samples (yk) are needed to 
achieve 0.99 probability of detection. 

Setting M = 100 (integration time = 1.25 sec), r 
is obtained using Eqs. (21) and (22), with Pj = 10 -4 . 
Figure 9 predicts that r should be set to 46, assuming 
that the outputs of the integrate- and- dump devices are 
scaled by the factor 1/ y/2cr%L. Using Fig. 7 (M = 30), 
one can check that when the detector is in-lock, SNR 2 & 
13 + 101og 1Q (100/30) = 16.3 dB, where scaling was per- 
formed to extend the results of Fig. 7 for M = 100. 
This is confirmed in Fig. 5, which depicts P d = 0.99 for 
SNR, « 16 dB. 


IV. Conclusion 

This article presents a mathematical model of the 
performance of two lock detectors for Costas loops: the 
square-law detector and the absolute-value detector. The 
model concentrates on the impact of phase jitter in the 
tracking loop on the performance of the lock detectors. 
Results of the analysis were verified by computer simula- 
tion and show that low loop SNRs result in a degradation 
in probability of lock detection, the amount of which is 
dependent on the scenario of interest. That decrease can 
be overcome by properly readjusting the design parame- 
ters. It was further shown that the square-law detector 
experiences less degradation due to phase jitter than the 
abolute- value detector, and that the degradation in de- 
tector signal-to-noise ratio is more pronounced for square- 
wave than for sine- wave signals. 
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LOCK DECISION 


Fig. 1. Costas loop with the square-law detector. 
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Appendix 

Derivation of Detector’s First Two Moments at Low Loop SNR 


A. Square-Law Detector 

Using Eq. (5) with Eq. (4), the expression for lock- 
detector signal is rewritten as follows 

Vk — L 2 Pd (w* - v\) + (n 2 Ik - riQk) 

-}- 2d k L\J~pD (wjbn/jt - v k riQ k ) 

= a-k + b k + c k 
M 

z = ^ (a* + b k + c k ) (A-l) 

Jt=l 

To assess the performance of the lock detector, the first two 
moments of z are needed. For a square-wave sub carrier, 
— v 2 — 1 — 7 | <j>k |) and for a sine- wave sub carrier, 
uj k — v k = cos 2<j ) k . Assuming that ^ is a zero-mean (no 
Doppler) Gaussian phase process, it can be shown that 


Eqs. (A-2) and (A-3) are independent of k . The expected 
value of z is 

/i, = M^L 2 P D d + La 2 n {l-W)) (A-4) 

where d is the signal degradation factor due to phase jitter 
in the tracking loop, and c* = 0. Using Eq. (A-2) results 
in 

d=i - 2 (S) ■ T ‘( i - exp (-^)) 

for a square-wave sub carrier, and 

d = exp(-2o-£) 

for a sine-w T ave subcarrier. To compute the variance of z, 
Eq. (A-l) is used to get 






' ( b 2 

^[i-expl-^y 


and 


<t>l = erf 


(A-2) 


l„+bexp(-±A (A-3) 


MM MM 

(»?.• - "«.) K - 4,0 

i= 1 ; = 1 i=l j = l 


M 


Af Af 


+ B ^ (wjtn/,- - rtriQi) 2 + C ^ ( n b ~ n Qj) 

1 = 1 j = l 


1 = 1 


where 6 = —2L, cr£ is the variance of the phase-jitter process 
in the subcarrier loop, and erf (x) is the error function 
defined as 

erf (ar) = -^= f exp (-t 2 )dt 
v n Jo 

At first glance, Eq. (A-3) seems to express the variance 
of the phase-error process _as a nonlinear function of itself. 
This is not the case since <j)\ is the variance in the window 
of the loop, and is the variance integrated over the 
complete density. This subtle effect is due to the Gaussian 
assumption, which approximates a density over a finite 
interval by the Gaussian density which is over an infinite 
interval. Numerically, <j>\ and are very close for all 
practical values of loop SNR. Note that the moments of 


where A — L 4 Pp> D — ALPPq, C — 2L 2 Po } and gi = 
(w? - v >?) = (1 - ~ | <f>i |) for a square- wave subcarrier. In 
the above equation, six terms were left out because their 
expected value is zero. Taking the expected value of z 2 , 
first over the thermal noise, and second over the phase 
process in the tracking loop, the following is obtained: 

F = L A PlM 2 g + 2ML 2 <t*{\ + W 2 ) 

+ M 2 l?cr A n (\ - W) 2 
+ AML 3 P D cr 2 (f + hW) 

+ 2M 2 L 3 al(l-W)d (A-5) 
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where 


M M 


= 772 


*=1 i=i 


and /, h are given by 


/ = w 


2 _ 

Jt — 



A = 





(A-6) 


for a square-wave subcarrier and 

f — u'* — cos 2 <j> = 0.5 ^1 + exp(— 2cr^)^ 

(A-7) 

h — sin 2 <f> = 0.5^1 — exp(— 2cr^)^ 

for a sine- wave sub carrier. The variance of 2 can be found 
from the relation cr 2 = I 2 - (z) 2 . Using Eqs. (A-4) and 
(A-5) results in 

<j 2 - M 2 L a P 2 d (g - d 2 ) + 4 ML 3 P D a 2 n (f + AW) 


to evaluate <7, one needs to know the second-order joint 
density function of the phase-error process <t>j , r )- As 
an approximation, it is assumed to be a twodimensional 
Gaussian density specified by the means, variances, and 
correlation coefficient 77 (t), which is obtained by simula- 
tion. Hence, 




27rv^ 2 (0) - R 2 {r) 

R( O)0 2 - 2fl(r)^j + R(0)<j>] 


exp -- 


2 (R 2 (0)- ft 2 (r)) 


- C 2 (r) 


exp -- 


<t> 2 _ - 2 C[r)Mj + 0) 

2<r 2 (l -C 2 (r)) 


Here g can be evaluated from the following: 


^ MM f f 

» = J j 9i9jP(^i,4>jy T ij) d< l , i d< l > i (A-ll) 


l 6 


(A-10) 


•'= 1 J = 1 -6 -i 


+ 2ML 2 <t*(1 + W 2 ) 


(A-8) 


r,j = T L {i - j) - ti - tj 


Note that at high loop SNR, g -* 1, d -+ 1, / -► 1, h -> 0, 
and Eqs. (A-4) and (A-8) reduce to 

Hz — ML 2 Pd + MLa 2 (l — W) 


a 2 = mi 2 a\ 


LP d 1 + W 2 


+ 


Expanding Eq. (A-ll) for the square-wave subcarrier re- 
sults in 

g= ^l-i|Tl+(i) ^ 

where 


as they should. In order to evaluate g, one must know the 
correlation between samples of the phase-error process in 
the tracking loop. That was obtained by simulation and is 
shown in Fig. A-l. A good closed-form model is given by 

R(t) = <rlC(r) (A-9a) 

where 

C(r) =(l- i||^) exp (-1.25B t r) (A-9b) 

and where <r£ is the closed-loop variance of the phase pro- 
cess and Bl is the one-sided loop bandwidth. In order 


-wiin I 4>i II 4>2 1 p{<i>\,^2,Tij)d4>id4>2 

(A-12) 


h b 


,=1 i =1 ~b -b 


The correlation function in Eq. (A-10) is symmetric, i.e, 
R(rij) = R(Tji), and depends only on the magnitude of 
the difference Jfc =| i - j |. This property allows the double 
summation in Eq. (A-12) to be reduced to a single sum- 
mation: 


Af-l 


= <*W*) 


k — 0 
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where c(0) = M , c(k ) — 2(M — k) for k ~ 1, 2, . . . , M — 1, 
and 



I tfi II <t>2 | P(^l,^2,rfc)^id<^2 


B. Absolute-Value Detector 

At low SNR in the tracking loop, the mean value of 
the detector’s signal z is obtained by taking the expected 
value of yk (Eq. 15) over the phase process in the tracking 
loop, and multiplying the result by M (the number of yk 
samples): 


Note that c (^) = M 2 as it should. For k = 0, the 

above probability density function (pdf) reduces to a delta 
function times a zero-mean Gaussian pdf with variance <rj, 
so that 

\ <t>i \\ h \ = 

Because 



the first term of Eq. (A-8) can be rewritten as follows: 

It can be shown that the lower bound of the above equation 
equals 

= ML ' P ° (;)’ (■ - ;) "5 

Keeping all detector parameters constant, as the loop SNR 
decreases (i.e., increases), fi z , cr 2 , SNR^, Pf } and Pd 
decrease. 

For a sine-wave sub carrier, 

1 Af-1 

t = JJJ £ 

k=0 

where c(k) is the same as in the square-wave subcarrier 
case, and 

{ 0.5 fl + exp(— 2<r^)^ for fc = 0 

exp(— 4 <t|) cosh ^2o-JC(rjt)j fc = 1, 2, . . . , M - 1 


/*2 = M(r Ik - r Qk ) 



where w and v are defined by Eq. (2) or Eq. (3), and 
r ik and rQk are defined by Eq. (14). The variance of z is 
obtained from ~~z 2 , namely 

MM MM 

I x n II X I] I + Y1 Y1 1 II I 

i = 1 j = 1 »=1 jzzl 

M M / M \ 2 

- 2^ Y! I Xl 3 II I - ( Y I X U I - I I ) 

i = l j = 1 \t = l / 

The following is now obtained: 

<j 7 z = M (l 2 P d + Lcr* (! + w)) 

+ Y1Y + 'WqJ) - M 2 i 7 ? + rQ 2 ) 

all i,j 

W 

where again rj and rg are defined by Eq. (14). Unfortu- 
nately, closed-form solutions for most of the above equa- 
tions are not obtainable and their evaluation has to be 
done numerically. 

The double sum in the above equation equals 
AT— 1 

s = ^ c(k)d(k) 

Jb=l 
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where c(k ) — 2 (M — k) for k — 1, 2, . . . , M — 1, and 



where the probability density function is again approxi- 
mated by Eq. (A-10). An upper bound of a 2 is 


M(L 2 P D + L*l( 1 + W) — (r7 2 +^ 2 )) 
and a lower bound is 

m(l 2 P d (/ + / 0 + + W)- (rj 2 + rg 2 )) 

where f and h are defined in Eqs. (A-6) and (A-7). 
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NORMALIZED CORRELATION 




